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ABSTRACT 

Numerical wave refraction programs permit a detailed study of the 
transformation of wave energy as waves move from deep water to shallow 
water. By eliminating the subjectivity that is present with hand drawn 
diagrams the effect of small variations in the initial assumptions and 
wave conditions can be investigated. The effects on wave refraction of 
tide level changes and bottom friction are investigated here. It is 
demonstrated that a uniform increase in water level as would be caused 
by tidal fluctuations can cause a significant change in the wave re- 
fraction pattern for a given nearshore region. A computing procedure 
is developed to permit numerical refraction programs to account for 
bottom friction. . The reduction in wave height caused by bottom friction 
depends primarily on bottom slope; this relation is shown in tabular 


and graphical form. 
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TABLE OF SYMBOLS 
approach angle of wave ray 


ray separation in shallow water divided by ray separation in 
deep water 


wave velocity 

deep water wave velocity 

dimensionless friction coefficient for the bottom. 
curvature of ray 

acceleration due to gravity 

wave height at position x 

deep water wave height 

wave height at position xj 

wave height at position x divided by deep water: wave height 


wave height at position x, when bottom friction is included, 
divided by deep water wave height 


wave height reduction factor due to bottom friction alone 
refraction factor 

combined friction and refraction factor 

direct shoaling factor 

wave length 

deep water wave length 

wave ray 

wave period 


distance measured along the wave ray in the direction of 
propogation of the waves 


x - XT 


3.1416... 


All values are assumed to be in the English System of Units. 


a Introduction 

Preparing accurate wave refraction diagrams is a necessary first 
step in many coastal engineering and amphibious operations. The advent 
of the numerical wave refraction programs [Griswold (1963) , Stouppe 
(1966) , Wilson (1966) | has considerably reduced the time required to 
construct a wave refraction diagram. The speed of the computer opera- 
tions permits more parameters to be included in the calculations and 
improves the accuracy of the output. Of equal importance is the fact 
that numerical products provide a completely objective basis for evalua- 
ting the cumulative effect of various assumptions on wave refraction. 

The present study has investigated the effects of tide-level 
changes and bottom friction on wave-refraction patterns for several 
environments. The investigation first compared the wave refraction 
programs of Griswold (1963), Stouppe (1966), and Wilson (1966). Since 
there was no way available to the author to compare program accuracy 
on an absolute scale, the comparison was based on method of computation, 
simplicity of use, and presentation of results. 

The program written by Stouppe was considered to have the most 
potential and was chosen for use in this investigation. The major 
advantages of Stouppe's program are these: 

(1) The program uses water depth directly as the interpolation 
surface for evaluation of wave speed, thereby limiting the 
number of computations and increasing the program versatility. 

(2) The program uses a second-order non-linear differential equation 
to determine the refraction coefficient at each point along 


the ray. Thus, the wave height at each point can be computed. 
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(3) The program uses constant time steps in the computation rather 
than constant distance steps. This reduces the distance be- 
tween computations as the shoreline is approached. Thus for 
the same number of computer operations a greater percentage of 
computations are in the shallowest area where refraction, 
friction, etc., are most important, and are changing rapidly. 

(4) The program includes a graphical output that plots the wave 
orthogonals, wavecrests and the shoreline. 

.. The following modifications were made to Stouppe's program: 

(1) The program was made to recycle, permitting more effective use 
of computer time when more than one diagram is to be constructed. 

(2) Provision was made to add a constant value to the depth field 

to investigate tide level effects. 

(3) A subroutine was added to account for bottom friction in the 


computation of wave height. 


12 


De: Computer Program. 
Depth grid. 

The first step in the utilization of this program is the construction 
of a grid of water-depth values for the desired area. The grid must be 
sufficiently large so that the starting point of all rays to be followed 
is in deep water, i.e., the ratio of water depth to deep water wave 
length (d/Lo) is greater than 0.5. By convention the x-axis is positive 
increasing toward the shore while the y-axis is positive to the left of 
the x-axis. The grid interval is selected so that the bottom contours 
are reasonably parallel to one another within a given square. 

In constructing the depth grid it is advantageous to record the 
average value of all depths within that grid square, rather than the 
value of water depth at the intersection point. This procedure tends to 
correct for small variations in the depth of water within the depth grid 
without necessitating the use of a smaller grid size. All actual depths 
are made positive; extrapolated depth values are continued on land for 
two grid units from the shoreline and are made negative. Beyond two 
grid units from the shoreline any arbitrary negative depth may be used. 
For depths on the shoreline itself zero is used. 

To reduce the time required to obtain the depth grid a clear plastic 
overlay with a black dot at each grid intersection point was prepared by 
initially locating the grid dots on white paper and producing a viewfoil 
transparency with an Ozalid copier, Model 400. 

Input data. 

Stouppe's order of data input was modified to permit the program 
to recycle. The present data deck consists of three sections. The first 
section contains only the value of MM and NN which designate the number 


of points in the depth grid in the x and y directions respectively. The 
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second section contains the water depth grid. These values are arranged 


so that the computer reads all the y values for each successive x-grid 


position. 


Water depths are read in fathoms and are converted to feet 


within the computer program. The last section contains one card for 


each set of input parameters for which it is required to construct a 


wave refraction diagram. The input parameters include: 


oe 


EGE 


TIDE: 


: HeaNT* 


Al: 


NOR: 


TIME: 


DIST: 


GRID: 


The grid position for the starting point of the first ray 
(feet). All starting rays must be more than two grid 
intervals in distance from the edge of the depth grid. 
The coefficient of bottom friction. 

The tide level, in feet, above chart depth for which com- 
putation is required. This value is subtracted from the 
depth grid after each refraction pattern is completed. 
The deep water wave height (feet). This is required for 
the friction subroutine. 

The initial wave period (seconds). 

The initial wave angle (degrees). The smaller angle 
measured between the positive Y direction and the initial 
wave crest. The angle is made negative if the slope of 
the wave crest is negative. 

The number of rays to be followed. 

The time interval between computations for advancing of 
the wave front (seconds). 

The distance between rays (feet). 


The grid interval (feet). 


This order of data input permits wave refraction diagrams to be 


computed for many combinations of input variables while requiring the 


computer. to assemble the program and read the depth field only once. 
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Computer output. 


The output from the computer consists of both a printed output and 
a graphical output. The printed output variables are these: 
XxX, The x and y coordinates for a given wave crest and ray 
number (yards). 
COREFR: The coefficient of refraction (K,) for the wave at the 
point x,™. 
HHO: The ratio of wave height to deep water wave height 
neglecting bottom friction. 
HHFR: The ratio of wave height to deep water wave height when 
bottom friction is included. 
NGO: Indicates that the ray has terminated (1), or it is 
continuing (2). 
DEPTH: The water depth for grid position X and Y (feet). 
The graphical plot of wave crests is programmed for the Calcomp 
160 system utilizing the DRAW subroutine programmed by J. R. Ward for 
the Fortran 60 system. See Appendix II, Stouppe (1966), for subroutine 
listing. The first and each succeeding third wave crest are plotted. 
The zero water depth points are used for contouring the shoreline on 
the graph. 
Appendix 2 contains a sample of the printed and graphical output 
for Monterey Bay, California. 
Computer operations. 
Since the details of the computer operations are discussed by 
Stouppe (1966), only the main features will be discussed here. The 


formulas used in the computations are listed in Appendix 3. 
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The water depth at the first point is computed by fitting the 
closest nine grid-point depths to a quadric surface by the least-squares 
method, An iterative procedure is then used to solve for wave velocity. 
It: is assumed here that wave velocity is a function of water depth and 
period only. 

The wave ray is then moved to the next point by solving for the ray 
curvature and projecting the ray forward a distance equal to the product 
of the wave speed for the point multiplied by TIME. At this new point 
the value of ray separation (4 ), coefficient of refraction Ca 
shoaling factor (Ks), H/Ho, and HHFR are calculated. 

This procedure is repeated until all orthogonals have been advanced 
one time interval. The new wave crest is plotted and the entire pro- 
cedure is repeated until all orthogonals have either gone off the grid 


or reached the” shoreline. 
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3. The Effects of Tide Level Changes on Wave Refraction Patterns. 


Differences in shoaling factor, wave speed, and ray curvature at specific 
depths. 


The effect on wave refraction of a uniform increase in water depth 
caused by tide level changes is cumulative and is dependent on the bottom 
depth contours. Therefore, to evaluate the effect properly, a wave ray 
must be followed along its entire distance from deep water to the shore. 
However, by calculating the value of shoaling factor (Ks), wave speed (C), 
and ray curvature (FK) at both high tide and low tide it is possible to 
evaluate the order of magnitude of the differences in the basic parameters 
of wave refraction caused by tide level fluctuations. 

The value of K,, C, and FK were calculated initially for 25, 50, 
and 100 foot water depths; second calculations were made for a water 
depth equal to the initial water depth plus a tide level increase of 
6, 10, and 16 feet. To evaluate these differences on a relative scale 


a percent difference was calculated by using the formula: 


Xr ~Xur 


(1) 
Xt 


| 
% difference X ; 


where 


= Varlabke represent ine uGamieemnsor PK 


Xpt = value of X at high tide water depth 


X1t = value of X at low tide water depth. 

The percent difference of K, was generally small and under certain 
conditions went to zero (for those values associated with the inflection 
point in the curve when K, is plotted against water depth). 


Figure 1 shows the percent difference in wave speed corresponding 


to the low tide water depth for a specific tide range and wave period. 
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It is apparent that the difference is often greater than 10 percent and 
increases with increasing wave period, increasing tide range and decreasing 
low tide water depth. 


The values of ray curvature (FK) were determined using the equation: 


PK > bi Aina BE ~ eh BE 7 @ 
where 
A = approach angle 
ac je = derivative of wave speed in the X and Y directions 
d 


respectively. 
Figure 2 shows the percent difference in FK versus tide range (computed 


using equation (2) with the following initial values: A = 30°, caro 
gC . AC for a twenty foot change in water depth starting with an 


OM 
initial depth of 100 feet). The values of percent difference plotted 


are representative of the percent difference for all values of 


af gt 
dH? OY 


the brackets in equation (2) have nearly the same value for high tide 


A, , and water depth; this occurs because the terms inside 
and low tide and tend to cancel out in the percent difference calculation. 


This reduces the percent difference calculation for FK to 


: Cig = Cue 


% difference FK = nee (3) 
Cny7 





The fact that the difference calculation is primarily dependent on tide 
range and wave speed is clearly shown in Figure 2. 

It is concluded that the change between high tide and low tide in a 
refraction pattern is primarily determined by the change in ray curva- 


ture (FK). 
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Difference in wave refraction patterns for specific locations. 


To investigate the combined effect of the changes in C, K, and FK 
caused by a uniform increase in water depth, wave refraction diagrams 
were constructed for several locations using both the low tide water 
depths and the high tide water depths. 

The first area considered was the California coast near Big Sur, 
between latitude 36° 14' 30" and latitude 36° 17' 42". This area was 
chosen because eight-second waves arriving on the coast will be refracted 
by the offshore bar (shown by the shaded area in Figures 3 and 4) at low 
tide but not at high tide... Thus tide level changes will cause selective 
refraction by the offshore bar. Figure 3 shows the refraction pattern 
for 8 second waves approaching from 270 degrees at low tide; Figure 4 
shows a different pattern for the same condition at high tide. | 

Figures 3 and 4 show two major differences in the refraction pattern. 
First, orthogonals 6 and 7 cross approximately halfway to shore at low 
tide but they do not cross at high tide.? Secondly, at low tide there 
is a convergence of wave orthogonals at Point A; this convergence is 


not present at high tide. 


Ito make the changes in the orthogonal pattern more prominent 
Figures 3-8 are plotted with an expanded ordinate scale. Also, the 
computer was not used to contour the wave crests but rather the wave 
orthogonals were contoured manually; the shoreline was indicated by a 
triangle plotted whenever water depth was zero. TIME equalled 30 
seconds for all diagrams. 


*Interpretation of the effect of crossed orthogonals is beyond the 


scope of this paper. In this discussion we will only be concerned with 
the difference in ray paths between low tide and high tide. 


Zu 


Table 1 shows the value of the wave-height ratio (H/H)) computed 
at four wave crests before the shoreline is reached for specific points 
along the Coast. This table clearly shows that the value of wave height 
can be considerably different at low tide and high tide. 

Figures 5 and 6 represent the same initial conditions as Figures 
3 and 4 except that the initial wave direction is from 310 degrees. 

The change in the orthogonal pattern is clear. It is noted that in 
Figures 3 and 4 the primary change in the orthogonal pattern (between 
low tide and high tide) is in the vicinity of points B and C, whereas 
in Figures 5 and 6 the greater change occurs near Point A. 

Wave refraction diagrams were also computed for the coast of Maine 
between Latitude 44° 39' 35"' and latitude 44° 44' 45". This area was 
chosen in order to determine if the combination of a large tide range 
(16 feet) and deep water depths (greater than 200 feet a short distance 
from the shoreline) would cause a change in the refraction pattern. 
Figures 7 and 8 show the refraction patterns for l6-second waves arriving 
from 070 degrees at low tide and high tide. The difference in refraction 
patterns is evident. Table 3 shows the H/Hp ratios (calculated at four 
wave crests before reaching the shoreline) for specific locations at 
low-tide and high-tide conditions. 

In total, about 20 sets of refraction patterns were analyzed and 
it appears impossible to predict what changes in wave height and re- 
fraction patterns will occur as the result of tide level changes. The 
only generalization that can be made is that the curvature of an or- 
thogonal is generally greater at low tide than at high tide. ft is 
concluded that precise water-level information, including tide state, 


should be incorporated into wave refraction procedures. 
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TABLE [ 
H/Hg RATIO. AT LOW TIDE AND HIGH TIDE FOR SPECIFIC COASTAL LOCATIONS 


NEAR POINT SUR, CALIFORNIA. 


LOCATION H/Ho RATIO 
(SEE FIG. 3) LOW TIDE HIGH TIDE 
A ia fe 2 
B 1.3 9 
G 3 al 1.7 
D .8 8 
TABLE II 


H/H, RATIO AT LOW TIDE AND HIGH TIDE FOR SPECIFIC COASTAL LOCATIONS 


NEAR EASTERN HEAD, MAINE. 


LOCATION H/Hg RATIO 
(SEE FIG. 7) LOW TIDE HIGH TIDE 
A 2.8 B43 
B 10 Law 
C af .8 
D | 5 
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4, Bottom Friction Effects on Wave Height. 
Introduction to bottom friction. 

The present method for Paeeiee ine wave height assumes that the 
propagation of energy between two adjacent orthogonals is conserved as 
the wave proceeds from deep water to shallow water. This implies that 
the addition of energy by surface winds, dissipation of energy by bottom 
friction or internal friction, and the reflection of energy by a steep 
‘slope are absent. In nature this is not the case - the wave eaetee 
between adjacent orthogonals may be changed by all of these means. 

According to Lamb (1932), in fluids of small viscosity (such as 
sea water), the velocity of the wave and the relationship between depth, 
length, and period are unaltered by internal friction. This is based on 
the assumption of a sinusoidal wave form but the result appears to in- 
dicate a very slow rate of damping of all waves in deep water, which is 
in accord with the Limited observations in nature (Putnam and Johnson, 
1949). 

At the sea bottom, however, friction can cause a loss of wave energy. 
It can be demonstrated both theoretically and experimentally that the 
particle velocities of shallow water waves are finite at the bottom and 
actually approach the theoretical orbital values at a short distance from 
the bottom. Therefore, a boundary layer in which mechanical energy is 
dissipated through the distortion of fluid elements must exist at the 
bottom. The Breck roughness of a sandy bottom probably is determined 
by the size of the ripples in the sand rather than by the size of the sand ~- 
grains. The ripple size, however, appears to be determined by the size 
and sorting coefficients of the sand in addition to the height and period 


of the wave and the depth of water. 


30 


Putnam and Johnson established the general magnitude for the dimen- 
sionless friction coefficient for the sea bottom to be 0.01. (This value 
was computed assuming initial conditions of 12 second waves, 4 feet high, 
traveling over a sandy beach with a ripple spacing of 5 inches.) Putnam 
and Johnson point out that since the friction factor only enters the 
equations for wave height to the first power, the friction factor 
could change by + 50 percent without appreciably changing the order of 
magnitude of wave-height reduction. 

Based on the following assumptions Putnam and Johnson were able to 
develop a procedure to calculate wave height when bottom friction is 
included. 

(1) The oscillatory motion at the bottom is sinusoidal and 

corresponds to the wave theory based on small amplitudes. 

(2) The friction coefficient does not vary with velocity and 

is approximately independent of water depth. This assumption 
amounts to choosing an average friction factor for the entire 
region over which the sea bottom affects the wave motion. 

(3) Perpendicular flow at the bottom resulting from percolation 

in the bottom material is negligibly small. 

(4) The sea bottom has a constant upward slope from offshore to 

the breaker point. 
Computing procedures. 

Bretschneider and Reid (1954) extended the original work of 
Putnam and Johnson and developed new equations to solve for wave 
height which included simplified equations for the special cases of 


constant water depth and constant bottom slope. 
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Bottom friction is of importance only on gently sloping beaches 
and in shallow water; thus it is possible, when using numerical wave 
refraction procedures, to calculate the reduction to wave height due 
to bottom friction by making the assumption that the water depth is 
constant between successive computations along the ray path from deep 
watér to the shoreline. " | 


’“*The wave height at each step is computed by the equation 


HH, Lt {0K ) eo Kea Ke PT w 


where 
Hea elemcelctraeed wave height, neglectane bottom rriction, Lor 
Ene point: 
Ke = the refraction factor between point n-1l and n. 


The coefficient of friction (Kg) is calculated by the equation (from 


Bretschneider and Reid) 


a (5) 
-~ AB ay + | 


7) ene ae K's is : 
Fg? NT tinh 2h 
K. = direct shoaling factor; 


A* 


K¢ is computed at each ray point with Ax equal to the distance 


incremented distance. 


along the wave ray between point n andn+1. It is assumed in the 


computations that the water depth between point n and n +1 is constant 


and equal to the depth at point n. 
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Results of bottom friction computations. 


The accuracy of the approximation that water depth is constant can 
be seen by a comparison of wave height obtained by this procedure to 
the wave height obtained by the more sophisticated and accurate pro- 
cedure used by Putnam and Johnson. This comparison is shown in Table 3 
(computed for wave period = 12 seconds, bottom slope = 5 feet). The 
comparison is very favorable along the entire distance of travel from 
the point that the wave first feels bottom to the breaking point. The 
difference is certainly within the order of magnitude of the error in- 
troduced by other assumptions. The difference in values is caused by 
the approximation to constant depth steps in the calculating procedure. 

Since this program uses constant time steps in moving each ray point 
toward the shoreline, the horizontal distance over which depth is assumed 
constant can be reduced by shortening the computing time interval (TIME). 
Table 4 shows the reduction (in percent) of wave height at the breaking 
point for different computing time intervals (with wave and bottom con- 
ditions the same as shown for Table 3 above). As expected, there is a 
slight increase in wave height reduction as the computing interval is 
decreased. It can be concluded that the value of TIME will affect the 
wave heights calculated by this procedure, but the effect is small. 

Realizing the limitations imposed by the initial assumptions of 
Putnam and Johnson and the calculating procedure employed here, the 
computer program was used to compute the wave height (with bottom friction 
included) for varying initial periods and bottom slopes. The results are 
shown in Table 5. (In arriving at these results it was assumed that the 
initial wave height was 5 feet and the friction coefficient was 0.01.) 


The reduction in wave height (in terms of percent of the wave height 
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TABLE III wiped: 
WAVE HEIGHT REDUCTION DUE TO BOTTOM FRICTION AS DETERMINED BY PUTNAM 


AND JOHNSON AND BY THE NUMERICAL WAVE REFRACTION PROGRAM. 


T¥= 12 sec. Sa f= .O1 TIME = 30 sec. 
Bottom slope = 1:300 Breaking depth = 11.5 Ft. 
DISTANCE FROM PERCENT REDUCTION WAVE HEIGHT 
WHERE h = Lo/2 FARRELL PUINAM & JOHNSON 
TO BREAKING POINT CONSTANT DEPTH STEPS - CONSTANT SLOPE 
PCT PCT PCT 
0.0 0%: 0 OM 
50.0 0.0 0.0 
62.0 Ore 0.2 
Vics 5, 0.4 0.4 
SON5 geal ra 
86.8: ae, 2a 
90.8 Sell oe 
Oar 0 4.2 4.5 
I580 6.0 oe 2 
oy Al oe 7 oe 
100.0 Orr 5 20.8 
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TABLE IV 


PERCENT REDUCTION IN WAVE HEIGHT FOR VARIOUS VALUES OF TIME 


TIME (SEC) WAVE HT. REDUCTION ( %) 
20 L950 
30 135 
40 L320 
50 17.0 
TABLE V 


WAVE HEIGHT REDUCTION FOR BOTTOM SLOPES OF 1:100 to 1:400 AND WAVE 


PERIODS OF 8 to 16 SECONDS (PERCENT) 


WAVE BOTTOM SLOPE 
PERIOD 1: 100 L200 1:300 1:400 
8 5.4 14.0 21085 28.6 
10 Dis2 iS .5 2037 27.0 
EZ 52 1320 Loe 24.8 
16 520 2aez 18.2 24.2 
AVERAGE D2 S72 Ie ole) Zong 
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without friction) is shown for the theoretical breaking depth (from 
H.O. 234). To obtain the percent reduction at the specific breaking 
depth, it was necessary to fit a smooth curve to the plotted points of 
wave height reduction at the three closest water depths and interpolate 
the wave height reduction at the breaking point. 

Figure 9 shows the average reduction (for all periods from Table 5) 
of wave height versus bottom slope. The curve is extended by a dashed 
line to the reduction amount found by Putnam and Johnson for a 1:10 
beach slope. 

The quantitative value of Figure 9 is controlled by the input 
accuracy of the coefficient of friction, the initial assumptions and 
the computation procedure. However, Figure 9 does show, qualitatively, 
the effect of bottom slope on wave height reduction due to bottom 
friction. Figure 9 is considered sufficiently accurate to be used in 
surf forecasting on sandy beaches to provide a first approximation of 


wave height reduction due to bottom friction. 
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Figure 9. EXPECTED WAVE HEIGHT REDUCTION AT BREAKING 


POINT DUE TO BOTTOM FRICTION FOR BOTTOM SLOPES OF 


1:10 to 1:400. 
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Br Conclusions. 

The numerical approach to wave refraction is a powerful tool in 
that it eliminates human subjectivity and permits the refraction com- 
putations to account for more factors than is possible with hand drawn 
techniques. 

A computer program that uses constant time intervals between com- 
putations has more suitability in accounting for the effects of friction, 
wind, reflection, etc., than a computer program that employs constant 
distance — 

Tide level changes can have a significant effect on wave refraction. 
This is especially true in areas where there is a large tide range or 
an offshore bar that can cause selective refraction. 

The changes shown by the tide level effects also point to the im- 
portance of having an accurate water depth grid. To avoid significant 
errors in the refraction pattern a great deal of care must be exercised 
in the preparation of the depth grid. 

Bottom friction can and does have an effect on the wave height. 

The numerical refraction program provides a convenient way to account 


for this reduction. 
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6. Recommendations. 

The course of this investigation has permitted a close look at wave 
refraction from many points of view. The numerical refraction program 
is capable of yielding wave refraction information with great speed and 
accuracy. (Computer time is approximately 5 minutes for a 20. ray diagram.) 
However, before the advantages of the computer product can be fully 
realized the author believes further investigations and improvements are 
needed in the areas of input data, computing procedures and field in- 
vestigations. 

The need for accurate depth information is essential to the compu- 
tation of accurate refraction patterns. It is anticipated that the aeoene 
efforts of the Naval Oceanographic Office and the Coast and Geodetic Survey 
in the area of color photography will improve the water depth information 
especially in areas where the bottom contours change considerably and 
in the shallow areas where it is difficult for hydrographic vessels to 
proceed, — 

The computer program can be expanded to account for more variables 
and to provide additional output information; in particular: 

(1) The program should account for the fact that in nature there 

usually exists a spectrum of wave periods and directions, 
rather than the single period and direction that is considered 
here. 

(2) For amphibious operations it would be advantageous if the 

program determined the breaking depth of the waves, printed out 
the height of the waves at the breaking point, and then contoured 


the breaking point on the output refraction diagram. 
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(3) Another desirable output for amphibious operations is a computer 
plot on the output diagram of the area where the wave height is 


‘aminimum, as this would be an ideal area to anchor’Llarge ships. 


7% + 
, 


(4) The various subroutines should be checked to ensure that each 
section is providing comparable accuracy. 
Field work is necessary to determine the value of the bottom friction 
- factor for various type bottoms. This is essential before accurate wave 
height. information can be obtained. Field investigations are also needed 
to determine the accuracy of the computer program; until now the numerical 
wave refraction procedures have been checked only by comparison with hand 


drawn results. 
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APPENDIX 1 
Computer Program for Wave Refraction using Fortran 63 
on the C.D.C. 1604 Computer 
esas Title: FRictT! 


Variables of Subroutine: 


HHQ....6.2.-...-Height of wave at point X, Y after bottom friction 
has been included. 


BHO! call oe Deep water wave length. 
-RLQ.....-2-.2.2--Wave Length at point X, Y. 
» CFF,.........-efriction factor between point n and n+l. 
_ ROFF........0.-(1 - Kp) 
TRCFF..........Total of (1-Kg¢) terms to point n. 


fiz 


i TCE Gaal 2. At Lotal friction»factor, tempointwn: 


evSummary of Subroutine: 


The subroutine calculates the value of Kr and wave height reduced 


by bottom friction (HHFR) for each point along wave ray. 


aGee Stouppe (1966) for discussion of remainder of program. 
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Sample Computer Input and Output 


Dib 6 
3506 
4706 
4006 
4506 
3006 
5156 
2156 
3706 
2506 
5506 
2706 
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5056 
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460 6 
3156 
4356 
3056 
5156 
4006 
3256 
2656 
560 ¢ 
4306 


APPENDIX 2 


Sample Input 


5106 
SD 6 
4156 
2156 
4456 
3156 
480-6 
3156 
B25 
2156 
550. 
32006 


4706 
3906 
3806 
3006 
3806 
3606 
445.6 
2156 
307. 
3056 
2406 
3106 


4256 
455 
3356 
2206 
3306 
420 « 
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2346 
2606 
4056 
5506 
2806 


355% 
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1406 
2106 
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1906 
4706 
6056 
2406 


3006 
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1106 
170-6 
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6206 
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2506 
8106 
5936 
636 
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FARRELL WAVE REFRACTION PROGRAM MONTER EY BAY 
THEME L* 90 ANELE = 4 clap. JD) eats: 


APPENDIX 3 


SUMMARY OF EQUATIONS 


Wave Speed (C) 
ee ae 
Cmmmeienh ( Ae) | 


Ray Curvature (FK) 
= 1 : 4) es CT 


Coefficient of Refraction (K,) 
The value of ray separation (B ) is calculated by solving the 


second order, non-linear differential equation: 


DB + pee a3 =) 





Ds# “Ds 

where 
= J. =o ee ee 
P= tmoaA ¢ Xo Aiw Be 3 


g + An Ag re - Zatwhevoh £ 3S + coo At FA a4e 
The above equation is solved by the finite difference method. This 
results in an equation for the Beta value at the n+l point in terms 
of the Beta values at the awe previous points. The equation to be 
solved is then: 

as (PD -2)A, + (4-240°)& 

3 Z+PD 
where 
D is the incremented distance 
P,q are as defined above 


ye} ie B » are the Beta values of the two previous points. The 


The coefficient of refraction is calculated by the relation: 


kK, =o Vo 


hy, 


Direct Shoaling Coefficient (K,): 


H,= 3.2519 - 12, 8150(%) + £8.82 ( Gr) 720. 725UE) 
+ N. 6 81S($)" . 


Wave Height Reduction Factor Due to Bottom Friction (ey: 


K. ! 





RP! Ay. +] : 
Ks 
where 
A’: = 
_ 64 7’ Ks 5 
D, ° 39° «5, ZA } 
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